INTRINSIC ENERGY IS A LOOP SCHUR FUNCTION 



THOMAS LAM AND PAVLO PYLYAVSKYY 

Abstract. We give an explicit subtraction- free formula for the energy function in tensor 
products of Kirillov-Reshetikhin crystals for symmetric powers of the standard represen- 
tation of Ulj{sln)- The energy function is shown to be the tropicalization of a stretched 
staircase shape loop Schur function. The latter were introduced by the authors in the 
study of total positivity in loop groups. 



1. Introduction 

The intrinsic energy function plays an important role in the path model for affine highest 
weight crystals [7]. The energy function is also related to charge statistic of Lascoux- 
Schiitzenberger on semistandard tableaux (see [Hj), which establishes a relation between 
one dimensional configuration sums arising in solvable lattice models and Kostka-Foulkes 
polynomials, cf . [21 El E] • 

Let B = Bi ■ ■ ■ <S) Bm be a tensor product of U'^isi^) Kirillov-Reshetikhin crystals, 
where each Bi is the crystal for a symmetric power of the standard representation. We 
identify Bi with the semistandard Young tableaux with row shape, filled with the numbers 
1, 2, . . . , n. Let b = bi ® ■ ■ ■ bm & B , and write xf for the number of r's in bi. The 
upper index {r — i + is to be considered as an element of Z/nZ. Our main result is the 
following formula for the intrinsic energy function Db of B. 

Let 6t = {t,t — 1, . . . ,1) denote the staircase shape of side-length t. 

Theorem 1.1. We have 

Ds{b) = mm I Yl ^I) \ ■ 

where the minimum is over all semistandard tableaux T of shape (n — l)5m-i, and entries 
in 1, 2, . . . , m. 

In the physical interpretation, each 6, represents a particle, and the intrinsic energy 
function Dsib) is defined as the sum of (^) local energies of interactions of particles. 
Theorem 11.11 thus has the following interpretation: each tableau T encodes a way for m 
particles to interact simultaneously, and intrinsic energy is equal to the minimum of these. 

In [H], motivated by the study of total positivity for loop groups, we introduced a 
generalization of the ring of symmetric functions, called loop symmetric functions and 
denoted LSym. In particular, we defined distinguished elements of LSym called loop 
Schur functions (see Section [2]). It is shown in [9j that the algebra homomorphisms from 
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LSym to M taking nonnegative values on (skew) loop Schur functions are in bijection with 
totally nonnegative elements of the formal loop group. 

Recall that the tropicalization of a subtraction-free polynomial /, is obtained by re- 
placing multiplication by addition, and replacing addition by taking minimums. Theorem 
II. H is equivalent to 

Theorem 1.2. The function Db is the tropicalization of the loop Schur function s^^^-^^^^^ ^ 
in the variables {x^^'*''}. 



Theorem 11.21 is an immediate consequence of Theorems 12.51 and 13.21 below. 

Theorems 11.11 and 11.21 are canonical in the sense that they correspond to the monomial 
expansion of a polynomial. That a piecewise-linear expression for Dsib) exists is already 
clear from the literature. However, the fact a subtraction-free formula exists (or equiva- 
lently the rational version 0^(6) of energy is a polynomial with positive coefficients) is not 
apparent from the definition of Dsib), even though the latter takes nonnegative values. 
(See also Remark [31 ) 



Example 1. Let n = 2 and m = 3. Then 
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corresponding to the following tableaux of shape 82 = (2, 1): 
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We use in our calculations a birational analogue of the combinatorial R-matrix. It was 
previously studied by Kirillov [6J in the context of the Robinson-Schensted algorithm, by 
Noumi-Yamada [121 IIBJ in the context of discrete Painleve systems, by Etingof [3] in the 
context Yang-Baxter equations, by Berenstein-Kazhdan [Tj in the context of geometric 
crystals, and by the authors [H] in the context of total positivity of loop groups. 

2. Loop symmetric functions 
Fix an integer n > 1 throughout. 

2.1. Loop Schur functions. Let (x-'^M be a rectangular array of variables. 

V / l<i<m, r&/nZ 

We recall from [9] the definition of the ring of loop symmetric function^ in the variables 
xf \ denoted LSym^. A detailed study of loop symmetric functions will appear in [TO] . 

For k > 1 and r G Z/nZ, define the loop elementary symmetric functions and loop 
complete homogenous symmetric functions by 



Cfc (Xi,X2,...,X, 
fc (Xl,X2 



, . . . , X^) 



E 

l<ji<i2<---<iA;<m 



Jr) (r+1) (r+k-1) 



h *2 ik 



^In [9], there are two such rings: the ring of whirl loop symmetric functions, and the ring of curl loop 
symmetric functions. We use the former here. Furthermore, we only use finitely many variables here. 
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By convention, el''' = /i^''' = for /c < 0, and e, 



(r) 



(r) 



1. Note that e 



(r) 



for k > m. 



We call the upper index the color. When all n colors are identified, that is x^-*'' = x^^* for 
all i and s,s' G Z/nZ, these functions specialise to the usual elementary and complete 
homogenous symmetric functions [15]. We define LSym^ to be the ring generated by the 



e^j^\ Although it is not immidiately obvious, the h^^^ lie in LSym^. In fact, both the e^^^ 
and the h^^^ are instances of distinguished elements of LSym^ called loop Schur functions. 

A square s = in the i-th row and j-th column has content c{s) = i — j. We 

caution that our notion of content is the negative of the usual one. Let p/z/ be a skew 
shape. Recall that a semistandard Young tableaux T with shape p/u is a filling of each 
square s & p/v with an integer T{s) G Z>o so that the rows are weakly-increasing, and 
columns are increasing. For r G 'L/n'L, the r- weight x^ of a tableaux T is given by 

We shall draw our shapes and tableaux in English notation: 
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For n = 3 the 0-weight of the above tableau is 
define the loop (skew) Schur function by 



{1)n2/^(1)x3^(2)^{2)^(2)^(3)^(3)/ (3)x2 



^1 



X 



T 



where the summation is over all semistandard Young tableaux of (skew) shape p/u. We 
have the following analog of the Jacobi-Trudi formula. 



Theorem 2.1. [9, Theorem 7.6] We have s^y, , = det(e 



Proposition 2.2. For any k > we have 



-0 "'k 



Ar-^)u{r~2) {r-2), (r-3) 
'1 "'k-1 "T ^2 "'k-2 



...±e 



ir-k),{r-k-l) 







0. 



Proof. Let pi denote the hook shape {k — i,l, 

{r—i—l) (^~*) 



1) of size k. Then ^^h^k-i 



Spk-i ' + ■5pfc_i+i is the sum of two loop Schur functions (one of which is zero if i 
OT i = k). This holds since any two terms in ef and h^^_ 







viewed as tableaux of 

column and row shapes, fit together to give a semistandard tableau of one of the two hook 
shapes, depending on the entry in their smallest boxes. It is clear that as we sum over all 
i, all the hook shape loop Schur functions cancel out. □ 



^M^(r-l) 



2.2. Staircases. For A; > 1 and r G 'L/n'L, define 
^i''^(xi,X2,...,x™) = ^ 

I={h<i2<---<ik} 

where the summation is over multisets / C {1,2,3, 
more than n — 1 times. Note that if A; > m{n — 1) we have 
t'^'' lies in LSym^, but we shall not need it for what follows. 



X 



(r-fc+l) 



, m} such that no number occurs 

It) 

— 0. It can be shown that 
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Lemma 2.3. We have 



oo 



i=0 \seI,/nZ seZ/nZ 

where the Ci in the above formula denotes the usual elementary symmetric function. 

Proof. Let x^^^-'x^-^^"'^'' • ■ . x^^'^"'^'^^-' be a term in h^i^\ Let J C J be the set of indexes which 
occur in / = {zi < ^2 ^ ■ ■ ■ ^ ik} more than n — 1 times. Then the coefficient of this term 
on the right is equal to X]a'cj(^1)''^'- This is equal to 1 if \ J\ = and to otherwise. □ 

Lemma 2.4. For each k, we have 

oo 

J2('^y'tH^i, . . . , x™)rt7-^) (x„ . . . , x^) = 0. 

i=0 

Proof. Follows immidiately from Proposition 12.21 and Lemma 12.31 □ 

Define a^''^(xi, xs, . . . , x„) = ^J^g xf ^xj''"^^ ■ ■ ■ x!''"'+^V^!l7*^(x2, . . . , x^). 
Theorem 2.5. Form > 2 and r G Z/nZ, we have 

(r) I \ 

_ (r) / N (r+1) / N ('r+m-2)/ n 

~ ^{n-\){m-\) • • • 5 ^mjO'(„,_i)(m_2) 1^2, . . . , X^j • • ■ Cr(„_i) (^X^-l, X^j- 

Remark 1. When all the colors are identified, that is, "x^f^ = -* for all i and s, s' G Z/nZ, 
Theorem 12.51 is a coarsening of a result of Jucis [5], see also \l5i, Ex. 7.30]. 

It is clear that Theorem 12.51 holds for m = 2, for then it states that s,[[l]^(xi, X2) = 
all\{xi,X2). We shall prove Theorem 12.51 in SectionlH 

3. Affine crystals 

3.1. i?-matrix. We shall use [Mj as our main reference for affine crystals. 

Recall that a Kirillov-Reshetikhin crystal of U'g{s\n) is the crystal graph corresponding 
to the highest weight module with highest weight proportional to one of the fundamental 
weights. An affine crystal is the tensor product of several Kirillov-Reshetikhin crystals. 
We shall restrict our attention to the set € of affine crystals that are tensor products of 
symmetric powers of the standard representation. Each element b E B of such a symmetric 
power can be identified with a single row semistandard tableau in the alphabet 1, . . . ,n. 

If Bi,B2 are Kirillov-Reshetikhin crystals, the combinatorial R-matrix is the unique 
isomorphism Rbi,B2 : -Bi ® -B2 — ^ -B2 ® -Bi of affine crystals. It is known that the combi- 
natorial i?- matrices generate an action of Sm on i?i ® • • • ® -Bm- 

The combinatorial i?-matrix has a convenient interpretation [H] in terms of semistan- 
dard tableaux and the jeu de taquin algorithm [15]. Let 61 ® 62 be an element of Bi® B2. 
Then 62(^1 ® ^2) = Ci ® C2 G -B2 ® -Bi where Ci, C2 are the unique pair of row shaped 
tableaux which jeu de taquin to the same tableau that bi and 62 jeu de taquin to, as 
follows: 

^ ^ since both jeu de taquin to 112 2 4 
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The action of the combinatorial i?-matrix can be exphcitly described as follows (see 
[1]). Let Rbi, 32(^1 ® ^2) = ci (g) C2 and let Xi \ \ si{x^i^), Si{x2^) be the number of 
boxes filled with r-s in bi,b2,Ci,C2 respectively, r = 1, . . . ,n. Then 

Sl{Xi^) = X2^ + 7«r+l(&l, &2) - /^rl&l, ^2) and Si{x2'') = x'f^ + Kr{bl, 62) - /«r+l(&l, ^2), 

where 



mm 

0<s<n- 



X2 



t=l 



and the indexes are taken in Z/nZ. 
Example 2. In the example above {x{^ ^'•^^ -^'■^^ 



n-l 
t=s+l 



•^1 



(1,2,0,1) 
1 + 2-1 



«;i(&i,&2) 
= 2. 



) — (1)0,1,0), (^2 \x^2 \x\^ \x^2 



min(l,2,3,3) = 1, ^2(^15 ^2) = min(2,3,3,3) = 2 and 51(0;^ 



3.2. Intrinsic energy function. In [7J an important function Hb,b' '■ B ® B' ^ Z, 
called local coenergy is defined for a tensor product of two affine crystals. If B and B' are 
Kirillov-Reshetikhin crystals, local coenergy has the following simple description in terms 
of tableaux [141 ■ Given an element b ^ b' in B ^ B' & €, form a two- row semistandard 
tableaux from b and b' as above. After that measure the maximal number of cells one 
can slide the top row to the left so that we still have a valid semistandard tableau. This 
maximal number of cells is the value of HB,B'{b, b'). 

Example 3. For 
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the local coenergy is 3, because 
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is semistandard while 



is not. 



-(t) _^ -(1+i) 



X 



H=s+1 -^1 



It is easy to see that H B,B'{b^b') = Ki{b,b'), since each term 
is exactly the number of cells for which the boxes with (s + l)-s in them would allow 
to slide. The local coenergy is known to remain the same under the action of i?-matrix: 

Hb',b = Hb,b' ° Rb',b- 

We now define the intrinsic energy function D b : -B — > Z ([?]) of an affine crystal S G (£, 
following [13j. For Kirillov-Reshetikhin crystals energy is zero. Let 6 = fei ® 62 ® • • • ® &m 
be an element of an m-fold tensor product B = Bi ® B2 ® . . . ® Bm G Ct. We define the 
intrinsic energy DB{b) to be 



B>Bib) = ^ HB„Bj{SiSi+i ■ ■ ■ Sj_2{bj-1, 



(1) 

l<i<j<m 

Although not obvious from this definition, intrinsic energy is preserved by the i?-action. 
Example 4. Let us compute the intrinsic energy of the following element 6 = 61 (g) 62 ® ^3- 



12 3 



12 2 4 



1 3 
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We measure H{hi,h2) = 1, iJ(&2,^3) = 2. We also have seen above the result of applying 
Rbi,B2 to the first two tensor factors, which allows us to measure if (51(62), ^s) = 2. 
Therefore ^^(6) = 1 + 2 + 2 = 5. 
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3.3. Product (summation) formula for intrinsic energy. In this section, we switch 
from piecewise linear functions to rational functions. The two worlds are connected via 
tropicalization: if / is a subtraction-free polynomial, we let trop(/) denote the tropical- 
ization of /, obtained by replacing addition by minimum, and multiplication by addition. 



We are given a rectangular array of variables x. 



(r) 



Jr) _ (r+l~i) 



(1) ,-.(2) 

) • • • 5 -^j 

. Define 



,m, r ^ Z/nZ, with columns 



bi = i^Y' , xl^' , . . . , ). It is very convenient to make the following change of variables: 



71—1 S 



n-1 



s=0 t=l t=s+l 

SO that trop(Kr.(6j, 6j+i)) = 'Er+j-iibj.bj^i). 

In the variables x^-'^'', the birational i?-matrix acts (see [ini Proposition 3.10) via algebra 
isomorphisms Si, S2, . . . , Sm-i of the field of rational functions in {Xj^^"*}, given by 



and sAx 



(r) 



Kr-.l{bj, bj+i) 



Kr{bj,bj+i) Kr{bj,bj+i) 

and Sj{x^i^^) = x^*"^ for k 7^ j,j + 1. We also have E{bj ® 6j+i) = Kj{bj,bj^i), the rational 
analogue of the local coenergy. We let denote the rational analogue of the intrinsic 
energy function, so that trop(DB) = Db- The main result of this section is a product 
formula for Dg. We remark that Kirillov [6j has also studied the rational functions H and 



Lemma 3.1. Suppose 1 < i < j < m. Then 



^r('SjSj+i ■ ■ ■ Sj_2(6j-l), bj) 



a 



(r-j+i) 
(n-l)(i-i) 



and 



SiSi+i ■ ■ ■ Sj-i[X 



■j ' 



(r-j+i) (r-j+i—l) / n {r-j+i) / 

X- (T^^_-,^^j^^._.^(^Xj, Xj+i, . . . , Xj) Crj.^_i)(,-_j)+i(,Xi, Xj+i, . . . , Xj 



() — j+i) / ^ 



(r—j+i) / ^ 



Proof. We prove the two statements in parallel by induction on j — i. For j — i = 1 they 
coincide with the formulae for the Kr and the i?-action of Sj. By the induction assumption 



(r—j+i+l) {r—j+i) / n 

^{n-l)(j-i-l)l^*' ^i+ly • • • ; Xj_iJ 



0" 



(n-l)(i-l-i) 



^Our variables x^'^'' are nearly the same as Yamada's , differing by a reversal of the orientation of 
the circle. 
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Therefore 

^r(SiSj+i ■ ■ ■ Sj_2(^i-l)) bj) 
n— 1 s 1 

s=0 t=l t=s+l 

n-1 U'^-j+i Ar+t) (r+s-j+i+l) , \ W 

ELLt=s-j+i+2 0(n-l){j-i-l)\-^^' • • • ' -^J-^J LLt=l -^j 

(r-j+i) 



(r—j-i-ij / \ 



s=0 

(r—j+i) / N 
_ Cr(n-l)0-j)l^*' • • • ' ^iJ 

(r—j+i) / \ ' 

The last equahty holds because in a term of a\^^J^^._,^-^{xi, Xj+i, . . . , x^) the number of the 

x^*^-s is at most n — 1, while the number of the xf^-s and the x^*''-s together should be at 
least n — 1. 

Now we can also prove the second claim, since 

/ (r)N SjSj+i • ■ ■ Sj_2(x^-_-^ ■')/tr_l(SjSj+i ■ ■ ■ Sj„2(&j-l)) ^i) 
SiSi+i ■ ■ ■ Sj_i(x} ') = - 



t^r{SiSiJ^l ■ ■ ■ Sj^2{bj-i), bj) 

(r-j+i) (r-j+i-l) / N a^'^'^^'^'^Ux- x) a^'"^^''^ (x- X- ,) 

■ ■ ■ ' "(n-l)(j-i)^ *' ■ ■ ■ '^3^ " (n-l)(j-i-l)y^^^ • • • ' ^3-^) 

J'-i+^) Z~) (r-i+i-1) 77 N (r-i+i) 77 77 

(r-i+j) (r-j+i-1) / X 
_ "(n-l)(j-i)y^^' ■ ■ - i^j) 



^[n-l)(j-i)i^ii • • • ' ^i) 



□ 



Theorem 3.2. We have 

^B{b) = 0-["^i)(„_i)(xi, . . . , Xm)cr[^^_^^(^_2)(x2, . . . , X^) ■ ■ ■ j^^^ (Xm_i , X^) . 

Proof. The result follows from Lemma [XT] and ([T]), using H^-^Bj (sjSj+i ■ ■ ■ Sj^2{bj-i)0bj) = 
i^j-i{U,j-i{bj.i),bj). □ 

Remark 2. Comparing [111 Theorem 4.1] with Lemma [3.11 and Theorem 13.21 one can see 
that the tropicalization of the cr|^^_-^-jj.^_-^_^^(xj+i, . . . , x^) is essentially the index ind(m— i) 
in the index decomposition of charge in [TT] . 

Remark 3. Theorem 13.21 gives the irreducible factorization of 0^(6) (whereas Theorem 
11.21 gives the monomial expansion). To see this, one first notes that o"|^^_-^^(xi, X2) has 

(r) 

a unique monomial which involves X2 , from which one deduces the irreducibility. Now 
suppose that cr|^^_-,^^|.^_-,^^(xi, X2, . . . , x^) factorizes non-trivially as the product fg. We 

may write / = axm ^'"'"^ + b and = c as polynomials in Xm ^'"'"^ '^'^\ where 
a, 6, c do not involve Xm"^""^'"'™"^'*'*. One verifies that none of the variables Xm di- 
vide (t|^^_j^^^^_-,^j(xi, X2, . . . , Xm) and every monomial which contains Xm is di- 
visible by the product Xm ^''''xm ■ ■ - xti ""^^^"^ Thus we have 
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f = a'xm ^^^"^ ^^""xm '^^ ^^■■■Xm ^^^"^ ^^^""^^ + 6, wheie o' is a polynomial Hot 

involving any Xm ■ It is easy to see that a' cannot be a unit. But we then have a non-trivial 



factorization a'c = a 



(n-i)(m-2)(^i' • • • ' ^rn~i), and wc may proceed by induction. 
4. Proof of Theorem 12.51 



We let Am denote the Jacobi-Trudi matrix for the dilated staircase Schur funtion 



-,('■) 

\n-l)S„ 



By adding extra columns of size to (n — 1)5^-1 we may assume that is 



a na X na matrix. (Specifically, a = \{n — l)(m — l)/n].) 
Example 5. For n = 3 we have 



A, 



(r) 
3 

^0 



,('■-1) 
'4 

,('■-1) 
-3 

,('■-1) 

'1 



^1 



V 



•^3 
(r) 

(r) 



<^3 
^^1 



64 

p(--2) 



-0 



'1 



Lemma 4.1. Suppose n < i < na. Then column i of Am is obtained from column i — n 
by shifting the non-zero entries down by n — 1. 

Let Bm denote the n{a + 1) — 1 x n{a + 1) matrix obtained by adding n extra columns 
to Am, so that Lemma [4. II is still true. 



Example 6. For n = 3 we have 

/ eg ^4 

.C'-i) .('■-2) 



54 



^('"-1) ^('■-2) ^('^-1) 

(r-l) (r-2) (r) (r-1) 
C]^ ^2 eg 



Jr-2) 



(r-2) (r) (r-l) 
"0 '^l ^^2 

p(r-l) 



p(--2) 
C4 

^^2 

J'-2) 
^1 



Let 



r 



_(r-l) 



'(n— l)m' '''"(n— l)m— 1 ' 

be a column vector with components in LSym. 



(r-2) 



(r-ri(a+l)+l) 



(n— l)m— n(a+l)+l > 



Proposition 4.2. T/ie vector Bm ■ T is the zero vector. 



Is) 



Proof. It follows immediately from Lemma 12. 4[ and the fact that 
k > {n — l)m. 

Let Bi^m be the square matrix obtained from Bm by removing the i-th column. It is 



as long as 

□ 



easy to see that det(-Bj, 



'(n-l)5„/(i) 



is a loop skew Schur function. 



Proposition 4.3. We have det(i?j 



^(r-i-l) 
(n— l)m— i+1 



det(A 
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Proof. Since det(-Bi_m) = ^^{n-i)5m/ii) iisver the zero polynomial in the e^^'''*, we deduce 

that considered as a matrix with ceofficients in the field Frac(Xj-*''), the matrix i?^ has 
maximal rank. There is thus, up to scaling, a unique solution to the equation ■v = 0. 
By expanding the determinant of the matrix obtained from by repeating a row, it 
is easy to see that v = (det(Bi,m), — det(-B2,m), • • • , ±det(-B„a+n,m)) is a solution. But by 
Proposition 14.21 so is T. Thus the two vectors are proportional, and it remains to check 
that the scaling coefficient is det{Am)- We have det(i?i,m) = Cm . . . Cm "''^^^ det{Am) , 
while at the same time T^^-Vj-m ~ ^rnT^^ ■ ■ ■ em~"'~^^\ and the statement follows. □ 

Proof of Theore'm \2.5[ We have already verified the case m = 2, so we suppose that 
m > 2, and by induction on m that 

for every s G Z/nZ. We calculate that 

(^•-1) / \ 

(n—l)m 

- \^ ,,(^-l),,C'-2) . . „(r-i)„(r~l) / ^ 

- ^1 ^1 ^1 ■^(n-l)<5™./(i)'^^2,X3, . . . ,X^+lj 

i=0 
{n—l)m 

= ^ xf^^^xf'^^ ■ ■ ■ xf "'-'det(Si,m)(x2, X3, . . . , Xm+l) 
i=0 

{n—l)m \ 

^{r-l) (r-2) _ {r-i) ir-i-l) . Y Y ) 

i=0 J 
= '^(n-l)m(Xl' ^2, • • • , Xm+l)-s[|^''_i)5^_^ (xs, X3, . . . , X„+i) 

where in the first equality we used the tableau definition of s'l^n-i)&m ' penultimate 
equality we used Proposition 14.31 □ 

Remark 4. It is clear from the theory developed in [9] that loop Schur functions are 
invariants of the action of the symmetric group via the birational i?-action. Thus we 
have demonstrated directly that the energy function is an invariant of this action. This 
property is not obvious from the definition we use. 
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